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ABSTRACT 

Injection of material from a core-collapse supernova into the solar system's 
already-formed disk is one proposed mechanism for producing the short-lived ra- 
dionuclides, such as 26 Al and 41 Ca, inferred from isotopic studies of meteorites to 
have existed in the solar nebula. This hypothesis has recently been challenged on 
the basis that the injection of enough supernova material to match the meteoritic 
abundances of 26 Al and 41 Ca would produce large, measureable, and unobserved 
collateral effects on oxygen isotopes. Here we calculate again the shifts in oxy- 
gen isotopes due to injection of supernova material in the solar nebula, using a 
variety of nucleosynthetic conditions of our own progenitor explosions. Unlike 
previous studies of this type, we also consider the effect of non-homogeneity in 
abundance distribution of the nucleosynthesis products after the explosion. We 
calculate the shifts in oxygen isotopes due to injection of sufficient supernova 
material to produce the meteoritic abundances of 26 Al and 41 Ca, and analyze the 
predicted shifts in detail for compatibility with meteoritic data. We find that the 
range in possible isotopic shifts is considerable and sensitive to parameters such 
as progenitor mass and anisotropy of the explosion; however, a small number of 
compatible scenarios do exist. Because of the wide range of outcomes and the 
sensitivity of isotopic yields to assumed conditions, it is difficult to constrain the 
supernova that may have led to injection of 26 Al in the solar nebula. Conversely, 
we argue that the existence of viable counterexamples demonstrates that it is 
premature to use oxygen isotopes to rule out injection of 26 A1 and 41 Ca into the 
solar nebula protoplanetary disk by a nearby supernova. 
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Introduction 



From isotopic studies of meteorites it is known that the solar nebula contained at least a 



dozen different short- l ived radionuclides, o r SLRs (see reviews by iMcKeegan fc Davisl 12003 



Meyer fc Zinnerll2006l ; IWadhwa et al.ll2007l ). Identification of the sources of these SLRs could 
greatly constrain the Sun's birth environment and processes acting during star formation. 
The half-lives of some of these isotopes are shorter than the timescales ~ 10 6 — 10 yr typically 
associated with star formation, so they must have been produced near the time and place of 
the Sun's formation. The SLRs 26 Al, 41 Ca and 60 Fe, in particular, cannot have been inherited 
from the Sun's molecular cloud in abun dances consiste nt with ongoing Galactic nucleosynthe- 
sis, and must have been late additions ( ]Jacobsenll2005l ). A leading candidate for the source of 
these and other SLRs is one or more core-c ollapse supernovae in the Sun's birth environment 



contaminating eithe r its molecular cloud (ICameron fc Truran 



Gounelle et al. 


2009 


), or i 


2010; 


Loonev et al. 


2006) 



1977 



2000 



Vanhala fc Boss 



Ouellette et al. 



2005 



2002 



2007 



20061 ). Another leading candidate is pr oduction of SLRs by irradiation 



(by solar cosmic rays, essentially), within the solar nebula (IGounelle et al.l 120011 ). Because 
this latter me chanism by itself is inadequate t o explain the abundance of 60 Fe in the early 
solar system ( Leva et al. 2003 ; Gounelle 20061). it is gene rally accepted that the source of 
60 Fe is core-collapse supernovae (e.g. IWadhwa et al.l 120071 ). although it is not clear whether 
th e source of 60 Fe is a single, nearby supernova, or many (possibly distant) supernovae (as 



in 



Gounelle et al.l 120091 ). 



The origins of the other SLRs are also debated. A correlation between 26 Al and 41 Ca 
has been observed in meteorites, demanding a common source for these two isotopes af- 
ter the formation of the s olar nebula, as in the so called "late-injection" hypothesis of 
Sahijpal fc Goswamil (119981 ) . It is not yet clear whether these two SLRs are correlated with 
60 Fe. If evidence for a corelation could be found, this would strongly suggest that 26 Al and 
41 Ca were injected by the same supernova or supernovae that injected 60 Fe. The lack of such 
evidence, though, leaves open the possibility that 26 Al and 41 Ca were created by irradiation 
within the solar nebula while 60 Fe was injected separately by one or more supernovae, into 
the Sun's molecular cloud or protoplanetary disk. 

The abunda nces of the SLRs alon e have not yet enabled a discrimination between these 
possibilities, but IGounelle fc Meiboml (120071 hereafter GM07) have proposed that the oxy- 
gen isotopic ratios of early solar system materials may be used to rule out certain hy- 
potheses. Specifically, they argue that if 26 Al and 41 Ca were injected by a nearby super- 
nova into the Sun's pr otoplanetary disk, suffici ent to produce the observed meteoritic ratio 
26 A1/ 27 A1 5 x 10 -5 ( jMacPherson et al.l 119951 ). then the oxygen isotopic ratio of the solar 
nebula would be considerably altered: solar nebula materials formed before the injection 
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would have oxygen isotopic ratios significantly different from later-formed materials. GM07 
calculated the shifts in oxygen isotopic ratios accompanying injection of supernova 26 Al into 
the Sun's protoplanetary disk, using the isotopic yields in bulk supernova ejecta calculated 



by iRauscher et all ( 120021 ). A robust prediction of the GM07 models is that 17 0/ 16 of 
pre-injection materials should be significantly higher, by several percent, than post-injection 
materials. Examples of pre-injection materials may exist in meteorites, or especially in the 



solar wind sample returned by the Genesis mission (IBurnett et al.ll2003l ). Since preliminary 



resul ts from Genesis suggest the Sun is not isotopically heavy in oxygen (IMcKeegan et al. 



20091 ). and because no such 1 O-rich (or 16 0-poor) components have been discovered in 



meteorites, GM07 rule out a supernova origin for the 26 Al and 41 Ca in meteorites. 

The purpose of this paper is to reproduce and refine the method pioneered by GM07, 
and to test the conclusion that 26 Al cannot have a supernova origin. GM07 originally consid- 
ered only bulk ejecta of spherically symmetric supernova explosions. We begin our analysis 
with this case, but make necessary refinements to the method, and use our current nucle- 
osynthesis models to predict the isotopic yields. We then expand on the analysis of GM07, 
calculating the isotopic yields by allowing the disk to intercept ejecta from different parts 
of the supernova explosion rather than a uniformly mixed total yield, and by examining 
anisotropic explosions. We also simultaneously consider the injection of 41 Ca into the disk. 

The paper is organized as follows. In §2, we outline the method used to calculate shifts in 
oxygen isotopic composition due to supernova injection of 26 Al and 41 Ca, including updates 
to the method of GM07. In §3 we describe the results of nucleosynthesis simulations we have 
carried out, to determine the isotopic yields in supernova ejecta under various explosion 
scenarios. We determine the inputs needed to compute the shifts in solar nebula oxygen 
isotopic composition. These shifts in oxygen isotope before and after injection are presented 
in 84, and in 55 we draw conclusions. 



2. Method 

2.1. Calculation of Isotopic Shifts 

The method of GM07 is fairly straightforward. They assume that meteoritic components 
that sample the solar nebula's starting composition, before the acquisition of 26 Al, can be 
identified and measured. Likewise, they assume samples after the acquisition of 26 A1 can 
be identified and measured. Any difference in the oxygen isotopic content between samples 
of those two groups would then constitute a shift in oxygen isotopes brought about by the 
injection of supernova material. The preditced shift in oxygen isotopes due to injection 
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of supernova material into the protoplanetary disk can then be compared to the actual 
difference in oxygen isotopes before and after. In practice, because the vast majority of 
meteoritic components sample the solar nebula after injection, GM07 assumed a "final" 
value for the solar nebula oxygen isotopes, and used the isotopic yields in supernova ejecta 
to predict the initial composition. Testing the supernova injection hypothesis thus amounts 
to finding meteoritic inclusions with this initial oxygen isotopic composition. Such inclusions 
should have evidence for no live 26 Al at the time of their formation, and should be among 
the oldest meteoritic inclusions. 

The earliest-formed solids in the solar system are widely accepted to be the calcium- 
rich, aluminum-rich inclusions (CAIs), both because they contain m i neral s that are the 

and because 
It 



first solids expected to condense in a cooling solar nebula f Grossman 
their Pb-Pb ages are the oldest measured, at 4568.6 Myr (IBouvier fc Wadhwal 120091 ) 



19721 



is worth noting that because many of the minerals in CAIs are condensates, their isotopic 
composition should reflect that of the solar nebula gas. The vast majority of CAIs have 
inferred ini tial ratios 26 A1/ 27 A1 5 x 10~ 5 or appear to have been isotopically reset at a 



later date (IMacPherson et al.lll995f ). Only in a handful of CAIs known as "FUN" CAIs 
(fractionation with unknown nuclear effects) has it been possible to set firm upper limits on 
the initial 26 A1/ 27 A1 ratio and show these CA Is did not contain live 26 Al when they formed 
(jFahey et al.l Il987t IMacPherson et al.lll995l ). Thus, CAIs overall reflect the composition 
of the solar nebula at an early time, and FUN CAIs possibly record the oxygen isotopic 
abundance before the solar nebula acquired 26 Al. 

To make more precise statements, it is necessary to quantify the oxygen isotopic compo- 
sition of the nebula and various components. The molar fraction of oxygen in gas and rock 
can vary, so the relevant quantities are the ratios of the stable oxygen isotopes, 17 0/ 16 and 
18 0/ 16 0. In the field of cosmochemistry, these ratios are commonly expressed as deviations 
from a standard, in this case Standard Mean Ocean Wate r (SMOW), which has 17 0/ 16 = 
3.8288 x 10~ 4 and 18 0/ 16 = 2.0052 x 10" 3 fb'Neillll986L The fractional deviations of the 
isotopic ratios from these standard values are <5 17 and 5 ls O, and are measured in parts 
per thousand, or "permil" (% ). [That is, <5 17 = 1000 x (( 17 0/ 16 0)/( 17 0/ 16 0) SM ow - !)■] 
It is also standard to report the quantity A 17 ~ <5 17 — 0.52 5 18 0, because this quantity 
is conserved during almost all che mical fractio nation processes. [More precisely, A 17 = 
ln(l + 5 17 0) - 0.52471n(l + <5 18 0) JmSCS).] 



It is clear that the final oxygen isotopic composition of the nebula, (<5 17 0, <5 18 0)', will 
depend on its starting composition (5 17 0, 5 18 O) , the composition of the supernova material, 
(S 17 0, 5 18 0)sn, and the mass of supernova material injected (relative to the mass of the disk). 
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It is straightforward to show that 

s 17 o> - 5 17 o = (5 17 o SN - 5 17 o ) , (i) 

1 + X 

where <5 17 0sn is the isotopic ratio of the supernova material injected into the disk, and 
x = M( 16 0)sN/^( 16 0)disk measures the mass of injected oxygen relative to the oxygen 
present in the disk (with a similar formula applying to 5 18 0'). In terms of the masses 
involved, 

_ M( 16 Q) SN M( 26 A1) SN M( 26 Al) disk M( 27 Al) disk 
X M(™A\) SN X M( 2 6Al) disk X M( 2 ?Al) disk X M( 16 0) disk ' 1 > 

Most of these terms are defineable. First, M( 26 Al) SN /M( 26 Al) disk = exp(+At/r), where At is 
the time delay between supernova injection and isotopic closure of the meteoritic materials, 
and r = 1.03 Myr is the mean lifetime of 26 A1. By definition, M( 26 Al) disk /M( 27 Al) disk = 
(26/27) x (5 x 10 -5 ), because sufficient 26 Al must be injected to yield the meteoritic ratio. 
Fina lly, the isotopic abundances in the solar nebula are known (the ratio 27 A1/ 16 is taken 



from iLodderd (120031 )). so we derive 



4.846 x 10~ 7 , a / \ 

X = [M^M)/M^O)] SN ex P( +A ^ ( 3 ) 

Note that x is independent of the mass of the disk, but it increases with At, since larger values 
of At imply that more supernova material had to be injected to yield the same 26 A1/ 27 A1 
ratio, thereby implying larger isotopic shifts in oxygen associated with this injection. 

Besides the time delay At, the major inputs needed to infer (S l7 0, 5 18 0)o are the isotopic 
composition (<5 17 0, <5 18 0)sn and ratio of 16 to 26 Al in the supernova ejecta, and the oxygen 
isotopic composition of the po st-injection sola r nebu la. GM07 used bulk abundances of 



supernova ejecta calculated by iRauscher et al.l (120021 ) for the first set of quantities. They 



also assumed that the oxygen isotopic ratios of the post-injection solar nebula matched the 
SMOW values of the present-day Earth: (5 17 0', 5 18 0') = (0% ,0% ). This assumption is 
the main reason why they concluded that the pre-injection solar nebula had to be 17 0-rich, 
as we now demonstrate. Rearranging equation [TJ yields 

5 17 O = 5 17 0' + x (5 17 0' - cFOsn) • (4) 

Supernova ejecta tend to be 16 0-rich; in the extreme limit, 5 17 0sn ~ — 1000%o- If S 17 0' ~ 
0%o also, then <5 17 O ~ +(1000x)%o- That is, <5 17 O is inferred to have been positive and 
pot entially quite la r ge if re > 10~ 2 . The isotopic yields of the supernova ejecta computed 



by IRauscher et al.l (120021 ) were consistent with such large values of x and 5 17 0sn < 0, 
leading GM07 to conclude that generally 5 17 0o > 0%o- Indeed, for progenitor masses 
15 — 25 Mq, GM07 inferred S 17 Oq ~ +35 to +220 %o- Since there are no early-formed 
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meteoritic components with <5 17 this high, and becaus e the oxygen isotopic c omposition of 
the Sun appears to be consistent with 5 17 ~ — 60 %o (IMcKeegan et al.ll2009l ). GM07 ruled 
out supernova injection of 26 Al and 41 Ca. This conclusion depends on a few key assumptions 
that we update below. We consider the starting composition of the solar nebula, and take 
into account the non-homogeneity of supernova ejecta. 



2.2. Solar Nebula Oxygen Isotopic Composition 



Oxygen isotopic ratios potentially can test or rule out the supernova injection hypoth- 
esis, but several caveats must be applied to the method of GM07. The first and most 
important correction involves the oxygen isotopic composition of the solar nebula immedi- 
ately before and after the injection of supernova material. GM07 assumed the post-injection 
composition was equal to SMOW; however, SMOW is widely understood not to reflect the 
oxygen isotopic ratios of the solar nebula immediately after injection. On a three-isotope 
diagram of <5 17 versus <5 18 0, the oxygen isotopes of planetary and meteoritic materials are 
arrayed along a mi xing line called the C arbonaceous Chondrite Anhydrous Mineral (CCAM) 
line discovered by IClayton et al.l (119731). After correc ting for isotopic fractionation by ther- 
mal and chemical processes, lYoung fc Russell! (119981 ) inferred a mixing line with slope 1.0 
in the three-isotope diagram, and so we will refer to this mixing line as the "slope-1" line. 
Today the oxygen isotopic composition of the Earth (SMOW) is widely recognized to reflect 
a mixture of an isotopicall y lighter rocky component (to which C AIs belong), and an isotopi- 
cally heavy reservoir (e.g. IClaytonl 120031 ; lYoung &: Russelllll998l ). It is very likely that this 



component is isotopically heavy w ater, with S 17 0,5 18 > +30 %o (IClayton fc Mayedalll984 



Lyons &; Young! 120051 ; iLyonsI 120091 1 . The existence of isotopically heavy water is supported 
by the discovery (in the primitive carbonaceous chondrite Acfer 094 ) of a poorly character - 
ized product of aqueous alteration, with S 17 ~ 5 18 ~ +180%o (jSakamoto et al.l 120071 ). 
Quite possibly this heavy water is the result of a mass-de pendent photodissociation of CO i n 
the outer solar nebula by an external ultraviolet source (ILyons fc Young||2005l ; ILyonsI 120091 ). 
The photodissociation can be isotopically selective because the different isotopologues of CO 
molecules can self-shield; C 17 and C 18 are optically thin and dissociate more completely, 
releasing 17 and 18 atoms that react with H2 to form isotopically heavy water, while 
the abundant molecule C 16 is more optically thick and does not as completely dissociate. 
The light CO molecule is eventually lost with the nebular gas. Whatever the source of the 
isotopically heavy component, SMOW only represents a late stage in nebular evolution, and 
does not represent the state of the nebula immediately after injection of supernova material. 



Applying the same reasoning, it is likely that the starting composition of the solar neb- 
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ula was lower (more O-enriched) on the slope-1 line than most CAIs. The majority ofCAIs 



tend to cluster near (<5 17 0, <5 18 0) « (-41 % , -40 %o), i.e., A 17 « -20.2% fsee IClavton 
20031 . and references therein); but many of the most primit ive and unaltered C AIs cluster 
near (5 17 Q, 18 Q) w ( -50%o, -50 % ), or A 17 « -24 % JScott fc Krotll200lh . Likewise, 
Makide et aD J2009ah report A 17 = -23.3 ± 1.9 %o for "mineralogically pristine" CAIs. 



CAIs also contain grains of hibonite, spinel, and corundum, which are among the first min 



erals expected to condense from a cooling gas of solar composition (jEbe 



fc Grossman 2000) 



and which are presumably even more primitive than CAIs themselves. iScott &: Krotl (120011 ) 
report that hibonite grains are also found to clus ter near (5 17 0,5 18 0) ~ (— 50 %o, — 50 %o), 
or A 17 = — 24%o, while iMakide et al.l f l2009bl ) observed 4 hibonite grains from Allende 
and Semarkona to have oxygen isotopes in the range A 17 = —32 to — 17 %o- They also 
found that spinel grains from the CV chondrite Allende had A 17 = —25 ± 5%o, and 
that corundum grains from the CM ch ondrite Semarkona clustered strongly in the range 
A 17 = -24 ± 2% . iKrot et al.l (boioh likewise report A 17 = -24 ± 2% for primitive 
CAIs and amoeboid olivine aggregates, which are also believed to have condensed from solar 
nebula gas. From these results we infer that (5 17 0', 5 18 0') ~ (— 50 %o, — 50 %o) in the solar 
nebula immediately after the injection of supernova material. 

Meteoritic and other samples also constrain the initial (pre-injection) oxygen isotopic 
composition of the solar nebula, and find it to be very similar. As described above, very firm 
and low upper limits to initial 26 A1/ 27 A1 exist for FUN CAIs that m ar k them as haying f ormed 
before the injection of 26 Al and 41 Ca flSahijpal fc Goswamilll998l ). IKrot et al.l (120081 ) have 
identified a fractionation line associated with the FUN CAIs with A 17 = — 24.1°/oo that 
passes through (5 17 0,<5 18 0) ~ (— 51 %o, — 52 %o)- Presumably the original isotopic composi- 
tion of t he nebula matched t hat of the Sun, which might therefore be measured by Genesis 



mission (IBurnett et al.l 120031 ) . Preliminary measurem ents can be interpreted as clustering 
on a fractionation line with A 17 ~= —26.5 ± 5.6%o ( McKeegan et al.ll2009l ). which would 



intersect the slope-1 line at (£ 1 7 0, £ 18 Q) ~ 
A 17 « -33 ± 8% [2a errors] ( TMcKeegan et al.ll201oh . 



5 6 %p, — 57 %p), and other analyses suggest 



From these results it seems likely that the original solar nebula oxygen isotopic com- 
position was near the Genesis preliminary result of (5 17 0, 5 ls O) ~ (— 54°/oo, — 53 %o), or 
possibly much lower along the slope-1 line. Subsequent reaction of rock with a 16 0-depleted 
reservoir then moved materi al along the sl ope-1 line to (5 17 0,5 18 0) ~ (— 41 %o, — 40 %o), 
where most CAIs are found (jClaytonl 12003k and references therein). FUN CAIs appear to 
represent an intermediate stage in this process, only partially evolved along the slope-1 line. 
To fix values, we will simply assume the solar system protoplanetary disk isotopic ratios 
started as (<5 17 0, <5 18 0) = (-60%o, -60 %o). 



-8 - 



The above discussion changes the criterion by which one can reject the supernova injec- 
tion hypothesis. Because GM07 assumed an initial solar nebula composition near SMOW, 
they concluded that pre-injection samples necessarily would have had S 17 > 0, and the 
lack of such samples in meteorites ruled out the hypothesis. But we assert that the su- 
pernova injection hypothesis can be ruled out only if injection of supernova material nec- 
essarily shift the oxygen isotopic composition of the solar nebula from a composition near 
(S 17 0, 5 18 0) ~ (— 60 %o, — 60 %o) to one far off the slope-1 line, or one on the slope- 1 line but 
with 5 17 > — 50 %o- In this way, the GM07 method of using oxygen isotopic constraints 
might still allow a test of the supernova injection hypothesis. 



2.3. Magnitude of Isotopic Shift 

There are at least three scenarios wherein the shift in oxygen isotopes following injection 
of supernova material can be consistent with the above constraints. From equation (TJ it is 
seen that even if the supernova ejecta and the protoplanetary disk differ in oxygen isotopic 
composition by hundreds of permil, the shift in oxygen isotopes may be small (< 1 permil) 
if the injected mass is small, so that x < 10 -2 . More precisely, if 26 Al and the other 
SLRs are injected by a supernova into the solar nebula disk, then the magnitude of the 
shift in oxygen isotopes will depend on the fraction of ejecta oxygen that accompanies Al. 
In this first scenario, O and Al may be significantly fra c tiona ted during delivery of the 



ejecta to the solar nebula. For example, lOuellette et al.l (120071 ) find that effectively only 



material condensed from the supernova ejecta into large (> 1 fim radius) grains can be 
injected directly into a protoplanetary disk. In the extreme event that the only grains that 
entered the protoplanetary disk were corundum (AI2O3) grains, the isotopic shifts in oxygen 
would be negligible (< 0.001 %o)- Or, if only 10% of the oxygen in the ejecta condensed 
into grains, and 90% remained in gas that was excluded from the disk, then the isotopic 
shifts in oxygen isotope (for a given amount of injected 26 Al) would be 10 times smaller 
than predicted by GM07. It is therefore not possible to determine the shifts in oxygen 
isotopes following injection into a disk without quantifying the degree to which O and Al 
are fractionated between gas and solids. In what follows, we assume no fractionation, as 
such a calculation is beyond the scope of the present investigation; but we consider dust 
condensation in supernova ejecta to be a very important effect, one that potentially could 
significantly reduce the predicted isotopic shifts. 

In the second scenario, the shifts in oxygen isotopes could also remain small if the in- 
jected material was simply higher than expected in 26 Al (or lower in O), so that again x < 



10 . The calculations of GM07 relied on the bulk abundances calculated by lRauscher et al. 
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( 120021 ). That is, GM07 assumed that the injected material uniformly sampled the entirety 
of the supernova ejecta. Such a uniform sampling is unlikely, as supernovae often do explode 
in a clumpy fashion and asymmetrically. It has long been understood that asymmetries or 
hydrodynamic instabilities may disrupt the stratifica tion of the progenito r star, but they do 



not result in large scale c ompositional mixing (e.g. IJoggerst et al.l 120081 ; iHungerford et al. 



2005a 



Frvxell et a 



nant (IHwang et al 



19911) . The X-ray elemental maps of the Cassiopeia A supernova rem- 



as thous ands of clumps o 



20041 1 dramatically demonstrate that massive stars are likely to explode 
mat erial, each sampling different burning zones within the pro- 



genitor. lOuellette et al.l (120 lOf ) have argued that this may be a near-universal feature of 



core-collapse supernovae; at the very least, observations do not rule out this possibility. So 
it is more than possible that the solar nebula received materials from only limited regions 
within the ejecta in which the 26 A1/ 16 ratio could have varied considerably from the average 
value for the ejecta. The non- uniformity of the 26 A1/ 16 ratio may be magnified if the star 
explodes asymmetrically, allowing explosive nucleosynthesis to proceed differently even in 
parcels of gas in the same burning zone. 

Finally, in the third scenario by which isotopic shifts may conform to measurements, 
x need not be small, and the isotopic shifts may approach 10 %o in magnitude, so long as 
the injection moved the composition up the slope-1 line by ~ 10 %o (i-e., the change in <5 17 
equalled the change in <5 18 0, both being < 10 %o), or down the slope-1 line by a comparable 
or even larger amount. A shift from an initial composition (S 17 0, 5 18 O) ~ (— 60 %o, — 60 %o), 
consistent with Genesis measurements of the Sun's composition, to (<5 17 0, 5 18 0)' ~ (— 50 %o, — 50 %o), 
consistent with primitive meteoritic components, would not conflict with the data. Alterna- 
tively, a shift from an initial composition (5 17 0, 5 18 O) « (-60 % , -60 % ), to (<5 17 0, 5 ls O)' « 
(-70%o,-70% ), or even (S 17 0,S 18 0)' « (-80% , -80%o), followed by mixing with the 
16 0-poor reservoir that moves solar nebula solids up the slope-1 line, would also conform to 
the data. 

In the next section we compute the isotopic yields in core-collapse super novae of various 



proge nitor masses, both in spherically symmetric explosions (as considered by lRauscher et al. 



20021 ) and asymmetric explosions. These calculations allow us to predict the oxygen isotopic 
composition (S 17 0, 5 18 0)sn and the ratio x of the supernova material at various locations 
within the explosion, to assess the range of possible isotopic shifts under the second and 
third scenarios. A supernova injection scenario would be ruled out, unless either the in- 
jection of material results in small overall shifts (i.e. the injected material contains a high 
26 Al abundance relative to oxygen, or vice versa, a low oxygen abundance relative to Al), 
or the oxygen isotopes are shifted along the 'slope-1 line', in which case shifts of up to ~10 
permil in either direction are allowed. 
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3. Isotope Production within Supernovae 



3.1. Numerical Methods 



We calculate the yields of oxygen isotopes, 41 Ca, and 26 Al in several core-collapse su- 
pernova scenarios, listed in Table [TJ These calculations explore four different progenitor 
models, with a range of explosion scenarios for each. We use a large set of thermally driven 
ID explosions with varying kinetic energies and delays, for a star of initial mass 23 M Q and 
a more restricted range of explosions for a 16 M , and 23 M Q , with the hydrogen envelope 
stripped in a case B binary scenario, and a single 40 M progenitor that ends its life as 
a type WC/O after extensive mass loss. We also examine a 3D explosio n of th e 23 M 
binary progenitor. Details of the simulations can be found in lYoung et al.l (120091 ) . The set 
of progenitor models we selected by no means samples the entire diversity of supernovae, 
but it represents a variety of cases across a large range of progenitor masses and explosion 
parameters. It is sufficiently diverse to make generalizations for behaviors that appear across 
all models. 



Progenitor models were produced with the TYCHO stellar evolution code ( lYoung fc Arnett 
!} 5j). To mode l collap se and explosion, we use a 1-dimensional Lagrangian code developed 



by lHerant et al.l (119941 ) to follow the collapse through core bounce. This code includes 3-flavor 
neutrino transport using a flux-limited diffusion calculation and a c oupled set of equations of 
state to model the wide range of densities in the collapse p hase fseelHerant et al.lll994t iFrver 



19991 . for details). It includes a 14-element nuclear network ( jBenz. Thielemann, fc Hillslll989l ) 
to follow the energy generation. To get a range of explosion energies, we opted to remove the 
neutron star and drive an explosion by injecting energy in the innermost 15 zones (roughly 
0.035 M ). The duration and magnitude of energy injection of these artificial explosions 
were altered to produce the different explosion energies. Our 3-dimensional simulation uses 
the output of the 1-dimensional explosion (23 M binary star, 23m-run5) when the shock 
has reached 10 9 cm. We then m ap the structure o f this explosion into our 3D Smooth Particle 
Hydrodynamics code SNSPH (IFryer et al.ll2006l ) by placing shells of part icles whose prop 



erties are determined by the structure of our 1-dimensional explosion (see iHungerford et al. 



2005bl ; IFryer et al.ll2006l for details). The 3D simulation uses 1 million SPH particles and is 
followed for 800 seconds after collapse. 

The initial intent of the 3D simulation was to cre ate a fully 3-dimen sional calculation 
of an explosion with a moderate bipolar asymmetry (lYoung et al.l 120061 ) . The interesting 
behavior of 26 Al in the explosion then prompted us to consider the composition of this 
material in relation to an isotopic enrichment scenario of the solar system. Both observational 
and theoretical evidence indicate that asymmetry is strong and ubiquitous in supernovae (e.g. 
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Table 1. Explosion Simulations 



Simulation 


progenitor 


Energy 


MR emn 


Delay 










10 51 erg 


(Mg) 


ms 


23e-0.8 


23 


M 


single 


0.8 


5.7 


20 


23e-1.2 


23 


M 


single 


1.2 


4.1 


20 


23e-1.5 


23 


M 


single 


1.5 


3.2 


20 


23e-0. 7-0.8 


23 


M 


single 


0.8 


3.2 


700 


23e-0.7-1.5 


23 


M 


single 


1.5 


2.3 


700 


16m-runl 


16 


M 


binary 


1.5 


2.06 


20 


16m-run2 


16 


M 


binary 


0.8 


2.43 


20 


16m-run4 


16 


M 


binary 


5.9 


1.53 


20 


23m-runl 


23 


M 


binary 


1.9 


3.57 


20 


23m-run2 


23 


M 


binary 


1.2 


4.03 


20 


23m-run5 


23 


M 


binary 


6.6 


1.73 


20 


40m-runl 


10 


M 


single 


10 


1.75 


20 


40m-run5 


40 


M 


single 


1.8 


4.51 


20 


40m-run9 


40 


M 


single 


2.4 


6.02 


20 
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Frver et all 120071 : lYoung fc Frverl 120071 : iHungerford et al.ll2005bl : lLopez et all 120091 ). In the 



situation of supernova injection the asymmetric model's primary utility lies not in modeling 
a specific event, but rather sampling a wide range of thermodynamic histories for material 
capable of producing 26 Al. Injected material is likely to sample only a small region of 
the supernova, meaning we can treat each SPH particle as an isolated trajectory whose 
further evolution is not dependent upon the progenitor star or or global parameters of the 
explosion. Any explosion that produces a similar thermodynamic trajectory will end up with 
similar yields. We can thus probe a large variety of explosion condition not accessible to ID 
calculations without a prohibitive investment of computational time. Therefore we consider 
this asymmetric simulation sufficiently generic to justify its usage for this investigation. We 
create an asymmetric explosion with a geometric aspect ratio and final kinetic energy axis 
ratio designed to be roughly consistent with the degree of asymmetry implied by supernova 
polarization measurements. To simulate an asymmetric explosion, we modify the velocities 
within each shell by increasing those of particles within 30° of the z-axis by a factor of 6; we 
will refer to these parts of the supernova as high velocity structures (HVSs). The velocities 
of the remaining particles were decreased by a factor of 1.2, roughly conserving the explosion 
energy. This results in a 2:1 morphology between the semimajor to semiminor axes ratio 
by the end of the simulation. We did not introduce any angular dependence in the thermal 
energy. At these early times in the explosion, much of the explosion energy remains in 
thermal energy, so the total asymmetry in the explosion is not as extreme as our velocity 
modifications suggest. For a detailed discussi o n on choosing asymmet ry parameters for 3D 
explosions see IHungerford et al.l ( 120031 . l2005al ) ; iFryer fc Warrenl ( 120041 ) . 

As noted above, although we are only considering one 3D model, the results we get 
from that model are representative of a range of nucleosynthetic conditions that may occur 
in multiple explosion/progenitor scenarios. Each parcel of gas follows its own density and 
temperature evolution, which is determined by the local velocity of the parcel of gas. It is the 
local conditions of the gas that matter; it is unaware of the global evolution. An example 
of the trajectories from the Ring and Bubble regions are shown in Figure [TJ The Figure 
shows the temperature, density, and radiation entropy evolution for representative 26 Al-rich 
particles in the explosion. The lines are labeled with 26 Al, 18 F, and 18 abundances at the 
emd of the simulation (before complete radioactive decay of 18 F). We see three classes of 
trajectories: high temperature and high entropy, high temperature and low entropy, and low 
temperature, high entropy. Predictably, the high temperature, low entropy trajectories tend 
to have low 18 F (and therefore 18 0) abundances due to the photodisintegration of 18 F into 
14 N + a. High temperature, high entropy trajectories have a higher reverse rate for that 
reaction, preserving slightly more 18 F. The low temperature particles have the highest 18 F 
abundance at the end of burning. 
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The asymmetric explosion samples trajectories with a large span of velocity evolu- 
ti ons reasonabl e for p lausible asymmetries. As we demonstrate with the ID models and 
in lYoung et al.l ( 120091 ). the sites of production for 26 Al are similar across a wide range of 
stellar masses, so as long as we sample the particle trajectories well in a single asymmetric 
explosion, the results are robust to very large changes in progenitor mass and explosion 
asymmetry. This assumption is valid, since we are looking for regions in the explosion that 
produce plausible abundances, not a bulk yield. 



The network in the explosion code terminates at Ni and cannot follow neutron excess, 
so to accurately calculate the yields from these models we turn t o a post-process step . 
Nucleosynthesis post-processing was performed with the Burn code (lYoung fc Fryerl 120071 ). 
using a 524 element network terminating at "Tc. The initial abundances in each SPH 
particle are the 177 nuclei in the initial stellar mode l. The network m achinery is identical 
to that in TYCHO (for details of the simulations see lYoung et al.ll2009l ). 



3.2. Production of 26 Al in ID and 3D explosions 

There are three primary sites for production of 26 A1 in a massive star and its accom- 
panying supernova. It can be produced by hydrogen burning at high temperatures in the 
shell-burning regions of massive stars or evolved AGB stars. But neithe r of these productio n 



sites is important to the supernova injection scenario, as discussed in lYoung et al.l ( 120091 ). 
In the ID simulations, the two dominant production sites are two peaks in 26 Al abundances 
that coincide with peak temperatures in the explosion of 2.2 x 10 9 K and 1.5 x 10 9 K, in 
material that has undergone hydrostatic C burning in the progenitor. The higher of the two 
temperatures is sufficient for explosive C and Ne burning, the lower of the temperatures is 
near explosive C burning. The production of 26 Al in both regions is due to a significant 
increase in the flux of free p, n, and a-particles. At higher temperatures, characteristic of O 
burning, 26 A1 is quickly destroyed. 

Within the 3D calculations, 26 Al is produced in two main regions (see Figure [2]), similar 
to the ID results. The first is a ring-like structure and an associated small bubble where the 
two HVSs emerge into a lower density region, and which we denote the "Ring" . Material in 
the Ring has undergone explosive Ne and C burning during the explosion at temperatures 
slightly above 2 x 10 9 K, and corresponds to the explosive C and Ne region in the ID sim- 
ulations identified above. The second region, further out at the terminal end of the HVSs, 
is denoted the "Bubble" (see Figure [2]). Within the Bubble, material has undergone hydro- 
static C burning and then experienced peak shock temperatures ~ 1.5 x 10 9 K during the 
explosion, and corresponds to the second of the 26 Al peaks in the ID simulations identified 
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above (which we will refer to as sub-explosive C burning region). While the production sites 
of 26 Al in the simulation in 3D occur in zones of about the same temperatures as in the 
ID cases, the peak in 26 Al abundance in those regions are reversed from the corresponding 
regions in ID (i.e. the peak that is higher in 26 Al in ID is lower in 26 Al in 3D, and vice 
versa). An important aspect of the Bubble is that due to the rapid expansion of the HVSs, 
its density drops rapidly, quenching some of the nuclear reactions. The decrease in density 
in the ID simulations occurred at a slower rate, conversely more of the 26 A1 was able to 
be processed into other species. This freezeout of nuclear reactions (suppressing subsequent 
destruction of 26 Al) in the 3D simulation is the reason for the higher production of 26 Al in 
the Bubble, as compared to the Ring. 

In 3D, the 41 Ca production occurs in only one main production site, in a region adjacent 
to and partly overlapping the Ring in the 3D simulation (see Figure [3]), and in both the 
explosive C/Ne and sub-explosive C burning regions in the ID calculations. The production 
of 41 Ca requires a high 40 Ca abundance as the seed nucleus, and the main production 
channel is p- and n-capture onto 40 Ca. The slower drop in density and temperature in the 
ID calculations tended to favor a low level production of 41 Ca, which is why its abundance is 
slightly higher as compared to the 3D calculation. In the 3D calculation, 41 Ca is produced 
in the Ring, but the faster expansion of the material there due to the velocity asymmetry 
shuts off the reactions faster than in the ID models, and the final 41 Ca abundance is lower 
than in ID. In the bubble region, the lower temperature and rapid density falloff preclude 
any significant 41 Ca production. 

Within the zones where 26 Al is produced, the M( 26 A1)/M( 16 0) ratios can differ signifi- 
cantly from the bulk abundances. For the ID models, these ratios can vary by a factor of ~ 1 
up to a factor of ~ 100 between the explosive C/Ne burning region and the sub-explosive C 
burning region, with a typical variation of a factor of ~ 2 — 3. For example, in model 23e-1.5 
the M( 26 A1)/M( 16 0) ratio varies from (6.7 — 8.4) x 10 -6 , and for model 16m-run2 varies from 
1.6 x 10" 7 to 1.7 x 10" 5 . These are to be compared to the abundances in the bulk of the 
ejecta, which are M( 26 A1)/M( 16 0) » 6.7 x 10~ 6 for model 23e-1.5, « 1.7 x 10" 6 for model 
16m-run2, and varies between 4.0 x 10~ 7 to 1.9 x 10~ 4 across all ID explosions. In the 3D 
model, the ratios are M( 26 A1)/M( 16 0) wl-4x 10~ 4 in the SPH particles in the Bubble, 
w 1 — 4 x 10 -5 in the SPH particles in the Ring, and 2.88 x 10 -5 for the bulk supernova 
abundances. Thus we see that injection of material from the Bubble brings in an order of 
magnitude less oxygen (essentially all 16 0) per 26 Al atom than injection of material from 
the supernova overall or from the Ring. 
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3.3. Production of O isotopes in 26 Al-producing regions 

The abundances and isotopic compositions of oxygen within a localized region of the 
supernova can vary significantly from the bulk values, as their production is sensitive to 
the density, temperature, and composition, and to their variati ons with time in that region. 



In the ID explosions, density falls off roughly as a power law (lArnettlll996l ). Because this 
maintains a high density in the region where 26 Al forms by explosive C and Ne burning, 18 O 
is effectively synthesized into heavier species. 17 is also synthesized into heavier species 
but is also created by neutron captures onto 16 0. The net effect is that both 17 and 18 
are reduced relative to 16 0, and the 18 0/ 17 ratio is reduced. In the 23e-1.5 model, their 
mass fractions in the 26 Al rich zones never exceed ~ 10~ 5 , and the isotopic composition in 
nearly all our ID cases approaches (— 1000 %o, — 1000 %o), effectively pure 16 0. 

The details of oxygen isotopic abundances in the 26 Al rich zones of the 3D explosion 
differ. The more rapid expansion of the material in the 3D calculation limits the processing 
of 17 and other isotopes into heavier species, so the yield of those is higher than in the 
ID calculations. As one of the main burning products of explosive Ne burning, 16 is quite 
abundant in the Ring. However, some of the free p and n produced during explosive burning 
capture onto 16 O, producing 17 0, so the Ring (as Figure [2] shows) is quite enriched in 17 
relative to the rest of the explosion. In the Bubble, 16 is not produced explosively, and 
is mostly left over from the progenitor. The increased flux of free particles also burns some 
of the 16 there to 17 and 18 0. The freezeout from the expansion limits the processing of 
these isotopes into heavier species, so the 3D explosion is richer in these isotopes than the 
ID simulation. 

Part of the reason for the large variation in 18 isotopic yields is that most of it is 
produced by decay of 18 F (ti/ 2 = 110 minutes), which was co-produced with 17 and 18 0. 
Thus it depends sensitively on how much 18 F is present once nuclear burning shuts off, 
which in turn depends sensitively on the trajectories taken by the gas. At low temperatures 
the classical decay reaction 18 F — > 18 + e + completely dominates, but at high temperatures 
(above ~ 10 9 K), and low proton de nsity, another deca y channel opens up for 18 F, and it 



can decay also via 18 F — » 14 N + a (IGorres et al.ll2000f ). The branching ratio of these two 
reactions is very sensitive to temperature at around 1 x 10 9 K, with higher temperatures 
overwhelmingly favoring the decay to 14 N + a. 

The amount of 18 F remaining at the end of burning is highly dependent on the time 
taken to drop below that temperature, and the density evolution, as high entropies favor 
the destruction over the synthesis. Because of the power law drop off, the density in the ID 
calculations stayed higher for a longer period of time, as compared to the 3D calculation, 
thus isotopes had a longer time window in which they could be processed to higher species. 
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The density of the 3D calculation dropped faster due to the increased velocities of particles 
to create the asymmetry, so the nuclear burning shut off earlier, and more isotopes like 17 0, 
18 0, or 18 F survived the nucleosynthesis of the explosion. This results in a substantial 
variation in the 18 F abundance between the ID and the 3D calculation, and that same effect 
(i.e. how quickly the density drops) is also responsible for the variation in abundance of 
particles in the Ring and the Bubble by the end of the 3D simulation. 

In the ID simulations the full decay of all 18 F after the explosion was calculated in the 
reaction network. The 3D simulation was terminated earlier in its evolution before complete 
decay of the 18 F. As the temperature at that point in the explosion was well below 10 9 K, 
we assumed that any 18 F still present would decay into 18 0, as this is the only significant 
channel at these lower temperatures. 



4. Solar system oxygen isotopic shifts accompanying 26 Al delivery 

4.1. Spherically Symmetric Supernova Explosions 

It is now possible to calculate the shifts in oxygen isotopic abundances before and 
after the injection of supernova material, using equation [TJ As described in §2, the initial 
composition of the solar nebula was probably close to (5 17 0, <5 18 0)o ~ (— 60%o, — 60 %o), 
and we adopt this as the starting value. Based on the numerical simulations of §3, we have 
calculated the ratios x and the isotopic abundances (<5 17 0, <5 18 0)sn within the ejecta overall, 
and within the regions where 26 Al is produced. 

We begin with the case of the ID explosions. The regions where 26 A1 is produced are 
those we identified as the C/Ne explosive burning region, and the sub-explosive C burning 
region. For our purposes, these regions were defined based on the 26 Al content. The exact 
amount of 26 Al produced varied among the simulations, but the (radial) abundance distri- 
bution of 26 A1 in each simulation showed two distinct peaks that were at least one order of 
magnitude higher than the average 26 Al mass fraction. Thus the 26 Al-rich regions in the 
ID simulations were defined to be at least one order of magnitude higher in mass fraction 
than the average distribution. The final isotopic composition of the solar nebula following 
injection of supernova material has been calculated first a ssuming; the material had th e av- 



erage (bulk) composition of the ejecta for comparison with iGounelle fc Meiboml ( 120071 ). and 
then that of one of these 26 Al-rich regions. The results are presented in Tables |2|- [4] and in 
Figure |H Our results for injection of bulk ejecta from ID explosi ons conforms closely to the 



findings of GM07 using the ID models of iRauscher et al.l (120021 ) . The ejecta are generally 



very 16 0-rich, with 5 17 and S 18 that are large and negative. 
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Table 2. Oxygen isotopic shifts followin; 


g injection from a 


ID supernova 






16m- 


23m- 


40m- 


23e- 


23e-0.7- 




average 


average 


average 


average 


average 








Bulk 






16Q 


1.97 M 


1.48 M 


3.29 M 


2.93 M 


5.44 M 


17 


2.24 x 10" 4 M 


7.67 x 10" 6 M 


1.80 x 10" 4 M 


3.6 x 10" 4 M 


2.94 x 10" 4 M 


18Q 


1.85 x 10^ 3 M 


1.53 x 10~ 5 M 


3.29 x 10" 5 M 


6.63 x 10~ 5 M 


1.27 x 10~ 4 M 


26 A1 


3.86 x 10^ 6 M Q 


1.91 x 10~ 4 M 


1.52 x 10" 5 M 


2.17 x 10~ 5 M 


2.15 x 10~ 5 M 


41 Ca 


4.06 x l(T 6 M Q 


2.63 x 10^ 6 M 


1.21 x 10" 5 M 


1.06 x 10~ 5 M 


1.49 x 10~ 4 M 


5 l7 


-7l9%o 


-987 %o 


-899 % 


-788 % 


-854 % 


5 18 


-58l% 


-996 % 


-997 % 


-993 % 


-990 % 


X 


0.247 


0.00375 


0.105 


0.0655 


0.123 


Final 5 17 


-19l%o 


-63.5 %o 


-139 % 


-105 %o 


-147% 


Final 5 18 


-163 %o 


-63.5 %o 


-149% 


-117%o 


-162 %o 


Final A 17 


-118%o 


-31.1 %o 


-65.6 %o 


-45.1 % 


-66.2 % 


At 


1.41 Myr 


0.65 Myr 


1.36 Myr 


1.31 Myr 


1.73 Myr 
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Table 3. Oxygen isotopic shifts following injection from a ID supernova 





16m- 


zom- 


4U111- 


AoG- 


OQo n 1 
zoe-u. (- 




q ^ 'ov q fro 
dvt;I dg,t; 


avei a^,e 


avei a^,e 


d "X rOV O (TO 
dV<3± d^C 


"x rni' fro 
dV<3± d^C 






explosive 


C/Ne burning 






16Q 


5.19 x 10" 1 M 


5.14 x 10- 1 M 


4.83 x 10- 1 M 


1.35 M 


1.86 M 


17 


1.47 x 10~ 7 M 


2.89 x 10~ 7 M 


2.76 x 10" 8 M 


8.79 x 10~ 7 M 


4.96 x 10" 6 M 


18Q 


4.41 x 10" 7 M 


4.02 x 10" 6 M 


1.70 x 10" 8 M 


1.64 x 10~ 8 M 


3.21 x 10~ 8 M 


26 A1 


8.98 x 1(T 7 M 


1.04 x 10~ 4 M 


4.55 x 10~ 7 M 


9.37 x 10~ 6 M 


1.53 x 10~ 5 M 


41 Ca 


1.74 x 10" 6 M 


2.74 x 10~ 6 M 


5.59 x 10^ 8 M 


6.70 x 10~ 6 M 


1.16 x 10" 4 M 


5 l7 


-999 o/oo 


-998 0/00 


-999 %o 


-997%o 


-9940/00 


5 ls O 


-1000% 


-996 0/00 


-1000%o 


-IOOOO/00 


-IOOOO/00 


X 


0.280 


0.00241 


0.514 


0.0698 


0.0587 


Final 6 17 


-266 o/oo 


-62.3%o 


-379%o 


-121%o 


-112%o 


Final 5 18 


-266 %o 


-62.2 0/00 


-379%o 


-121%o 


-112%o 


Final A 17 


-1470/oQ 


-30.6 0/00 


-2270/00 


-61.3 0/00 


-56.1 0/00 


At 


1.51 Myr 


0.77 Myr 


1.03 Myr 


1.34 Myr 


1.75 Myr 
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Table 4. Oxygen isotopic shifts following injection from a ID supernova 





16m- 


23m- 


40m- 


zoe- 


zoe-U. 1- 




average 


average 


average 


average 


average 






sub-explosive C burning 






16Q 


1.78 x 1IT 1 M 


3.85 x lO^ 1 M 


2.47 M 


1.84 M 


1.91 M 


17 


1.81 x 10" 4 M 


1.03 x 10" 6 M 


1.88 x 10~ 7 M 


4.18 x 10" 5 M 


4.86 x 10" 5 M 


18Q 


3.12 x 10~ 4 M 


7.32 x 10^ 6 M 


6.51 x 10^ 8 M 


8.31 x 10~ 5 M 


9.14 x 10~ 5 M 


26 A1 


2.95 x 1(T 6 M 


6.77 x 10~ 5 M 


1.54 x 10^ 5 M 


1.42 x 10~ 5 M 


8.75 x 10~ 6 M 


41 Ca 


1.65 x 10^ 7 M 


1.35 x 10^ 6 M 


8.42 x 10^ 6 M 


7.78 x 10~ 6 M 


8.44 x 10~ 6 M 


5 l7 


1503 % 


-990 % 


-1000% 


-944 % 


-938 %o 


5 18 


-221 % 


-995 % 


-1000% 


-980 % 


-979 % 


X 


0.0292 


0.00276 


0.0775 


0.0628 


0.106 


Final 5 17 


-15.7% 


-62.6 %o 


-128 %o 


-112 %o 


-144%o 


Final 5 18 


-64.6 %o 


-62.6 %o 


-128 %o 


-114%o 


-148%o 


Final A 17 


+19.2 % 


-30.7% 


-64.9 % 


-55.3 % 


-71.4%o 


At 


0.90 Myr 


0.72 Myr 


1.29 Myr 


1.29 Myr 


1.39 Myr 
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GM07 likewise found, using the ID models of iRauscher et al.l (120021 ). that the ejecta 
are depleted in 17 0, and in most cases 18 as well. (Their 15 M & case is enriched in 18 0, in 
contrast to our 16 ca s e). Si milarly to our 23 M Q models, both the 21 M Q and the 25 M Q 
model of IRauscher et al.l (120021 ) are very depleted in 17 and 18 0, and the composi t ion o f 
those ejecta approach — 1000 %o for both 5 17 and 5 18 0, although the IRauscher et al.l (120021 ) 
models tend to be sli ghtly richer in 17 and 18 than ours. The similarity between the 
Rauscher et al.l (120021 ) bulk abundances and those calculated for this work are unsurprising. 
Post-He burning stages rapidly destroy 17 and 18 0, resulting in the oxygen in interior zones 
being nearly pure 16 0. A 20-25 M^model has a very larg e oxygen mantle. The slightly higher 
fraction of heavy isotopes in the IRauscher et al.l (120021 ) arise in a somewhat larger He shell 
that results from a less accurate treatment of mixing in their earlier stellar models. 

These differences reflect the variability inherent in calculations of nucleosynthesis in 
massive stars, especially where small shifts in stable isotopes are concerned. We also find, as 
did GM07, that the isotopic shifts associated with injection from ejecta from ID supernova 
explosions tend to be large (tens of permil), but not in all cases. In the 23m runs, 26 A1 is 
produced more abundantly, and the 26 A1/ 16 ratios yield x < 10~ 2 in these explosions. In 
the 23m cases, the ejecta are particularly 16 0-rich, but relatively less oxygen needs to be 
injected per 26 Al because more 26 Al is produced. (It should be remembered that 23m is 
a binary case, where a significant fraction of the 17 and 18 have been removed by mass 
loss from the He shell, and production of 26 Al has been enhanced by higher peak shock 
temperatures relative to the 23 M single star models.) The isotopic shifts associated with 
injection from 23m ID explosions are typically < 3%o- 

For all the cases the shifts in both <5 17 and <5 18 are negative and similar in magnitude; 
that is, the injection of material from this supernova moves the composition down the slope-1 
line. Later nebula evolution would produce materials that move back up this line. The most 
favorable case is the 23m supernova. The magnitude and direction of the isotopic shifts 
associated with the 23m case are such that current measurements of solar nebula materials 
do not rule out this possibility, even for bulk abundances. As we argue below, however, bulk 
abundances are not the best representation of the abundances of injected material. 

Tables 121 El and S] and Figure S] also show the isotopic shifts associated with injec- 
tion from only 26 Al-rich regions within the supernova. Even when considering injection of 
26 Al-rich regions only, the conclusions are not much changed: the isotopic shifts in oxygen 
generally are many tens of permil, and make the solar nebula more 16 0-rich. The sub- 
explosive C burning region of the 16m model is the only case that does not move the nebular 
composition along the slope 1 line. In general the sub-explosive C burning in the higher 
mass models provide the best results, as they produce the highest ratio of 26 Al to oxygen. 
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Tables El [3J and H] also give the yields of 41 Ca produced in each of the ID explosion 
scenarios. The post-injection 41 Ca/ 40 Ca ratio is generally more than sufficient to match the 
meteoritic ratio, and a time delay is implied before isotopic closure, so that 41 Ca can decay. 
For the ejecta from the 23 M progenitors, the implied time delay (for 41 Ca to decay to a 
level 41 Ca/ 40 Ca = 1.4 x 10~ 8 ) for the 23m cases is ~ 0.7 Myr for all three regions considered 
(i.e. bulk, explosive C/Ne burning, and sub-explosive C burning). The implied time delay 
for injection from the 23e cases is ~ 1.3 Myr for all three regions, and from the 23e-0.7 cases 
is 1.4 Myr (for the sub-explosive C burning region) -1.7 Myr (for the other two regions). 
The time delays for the other progenitor cases are all within those range s. The range of these 
time d elays are very similar to the range of 1.0 - 1.8 Myr calculated by lGounelle fc Meibom 
( 120071 ) . The effect of this time delay is to cause 26 Al to decay, too, before isotopic closure, 
and to increase the isotopic shifts in oxygen. The shifts are increased by factors of 2 (for the 
23m) at the low end to 5.2 (for the 23e-0.7 cases) at the high end. If isotopic closure is to be 
achieved in a few x(10 5 - 10 6 ) yr flMacPherson et al.lll995l : Kita et all 120051 ). then injection 
from the 23e and 23e-0.7 would seem to introduce too much 41 Ca to match constraints. 
Injection from the more energetic 23m progenitor cases are consistent with a small shift in 
oxygen isotopes downward along the slope-1 line, as well as the final 41 Ca/ 40 Ca ratio of the 
solar nebula. 



4.2. Asymmetric Supernova Explosions 

In Table [5] and Figure HJ we present the yields of 26 Al, 41 Ca, and oxygen isotopes 
in various regions of the ejecta in our simulation of the 3D explosion. We calculate the 
isotopic shifts if the injection uniformly samples all of the ejecta (bulk), if it samples the 
Ring material, and if it samples the Bubble material. By design, membership in the Bubble 
and Ring material is defined by high 26 Al content. These 26 Al- rich regions do not have 
well-defined edges, instead fading out monotonically in 26 Al-abundance as one moves out 
into the surrounding ejecta (see Figure [2]). In order to not impose an arbitrary geometry on 
these regions we determined membership by 26 Al amount per SNSPH particle. We used two 
different lower limits or thresholds for inclusion - 1.5 x 10 -13 M Q of 26 A1 per particle for a 
maximum extent of the 26 Al rich region, and 1.5 x 10" 11 M per particle for a minimum 
extent of the 26 Al-rich region ("high 26 Al" case). The most 26 Al-rich SPH particles in the 
Ring and Bubble had 1.5 x 10~ 10 M and 4.8 x 10~ 10 M of 26 A1, respectively. Each threshold 
picked out all SPH particles in the respective regions that it identified. 

Overall, the ejecta of the 3D simulation are much richer in 17 and 18 than the 
ID simulations, but also contain two regions (the Ring and the Bubble) in which the 26 Al 
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Table 5. 26 Al and O in the 3D explosion 





Bulk 


Ring 


Ring 


Bubble 


Bubble 








high A126 




high A126 


26 Al ( M ) 


1.474 x 1(T 6 


3.510 x 10~ 7 


2.420 x 10~ 7 


9.996 x 10" 7 


9.585 x 10~ 7 


16 O(M ) 


0.511 


2.888 x 10~ 2 


5.949 x 10~ 3 


7.025 x 10" 3 


2.550 x 10~ 3 


17 O(M ) 


2.170 x 10~ 4 


3.400 x 10~ 5 


2.241 x 10~ 5 


1.021 x 10" 5 


1.738 x 10~ 9 


18 O(M ) 


1.509 x 10- 1 


3.553 x 10- 3 


1.018 x 10" 4 


7.314 x 10~ 3 


5.635 x 10- 3 


41 Ca(M ) 


2.132 x 10~ 8 


2.039 x 10~ 8 


4.312 x 10~ 9 


1.827 x lO" 11 


1.239 x lO" 11 


5 17 


+43.9 % 


+ 1894% 


+8258 % 


+2573 % 


-998.3 % 


5 18 


+ 129887 % 


+53546 % 


+6582 % 


+460545 % 


+978646 % 


X 


0.16808 


0.03987 


0.01192 


0.003406 


.001289 


Final <5 17 


-45.1 % 


+14.8 %o 


+37.9 % 


-51.1 %o 


-61.2 % 


Final 5 18 


+ 18638 % 


+ 1995 % 


+18.2 % 


+1503 % 


+ 1200 % 


Final A 17 


-1608 %o 


-561 % 


+27.8 % 


-534% 


-477 % 


At 


0.66 Myr 


0.90 Myr 


0.70 Myr 
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production is increased over the ID calculations. As we have previously discussed, the 
production of 18 is significantly altered from the ID results. The added yield from the 
decay of 18 F to 18 makes the ejecta significantly richer in this isotope, and results in large 
(tens of permil) to very large (hundreds of permil) and positive shifts in 5 18 for material 
from both the Ring and Bubble, and the bulk. This is in stark contrast to the 18 poor ejecta 
produced in the ID simulations, and emphasizes the sensitive dependence on the prevailing 
thermodynamic conditions of 18 O production. The production of 17 is much less sensitive 
to the thermodynamic conditions. In the 3D simulation we also see an increase over the ID 
cases in the production of 17 in the 26 Al rich regions and the bulk; however the change is 
not as drastic as in 18 0. This again differs from the ID calculations, and the more 17 0-rich 
ejecta result in positive shifts in 5 17 0, on the order of —1 to + 15 %o for the Bubble and bulk, 
and close to +100 %o for the Ring. 

Table [5] also shows the Ca/ 40 Ca ratio following injection of material from the 3D 
supernova into the solar nebula. If injection comes from the Bubble region only, the amount 
of 41 Ca injected is too low to conform to meteoritic ratios, and injection from the Bubble 
can be ruled out on these grounds. Injection of material from the Ring or bulk regions, 
in contrast, imply reasonable time delays ~ 0.66 — 0.90 Myr. This implies an increase in 
oxygen isotopic shifts of < 2.5 over what is presented in Table [51 These time delays are just 



below the ones iGounelle Sz Meiboml (120071 ) calculate, which again is explained by the faster 



density- drop off in the 3D calculation producing slightly less 41 Ca than in ID. 

When an explosion samples a variety of thermodynamic trajectories through asymmetry, 
including those that result in freeze-out conditions due to rapid expansion, the overriding 
conclusion to be derived is that a very large range in oxygen isotopic shifts is allowed. It 
would seem extremely unlikely that conditions in an asymmetric explosion would conspire to 
yield a small isotopic shift consistent with the meteoritic constraints, though more "normal" 
trajectories that do not experience this freeze-out process are still candidate production sites, 
as wee see in ID. 



5. Discussion 



As iNichols et al.l (119991 ) strongly advocated, injection of supernova ejecta can produce 
measurable "collateral damage" to stable isotope systems in protoplanetary disks. GM07 in 
particular point to the role of oxygen isotopes in constraining this process. The point of that 
paper was that the injection of 26 Al (and 41 Ca) from a single nearby supernova necessarily 
would have brought in significant levels of oxygen isotopically distinct from the pre-injection 
solar nebula. The solar nebula after injection, they argued, would differ in its oxygen isotopes 
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by several tens of permil from the pre-injection values, which they robustly predicted would 
be more 17 0-rich than the solar nebula. They cited the Genesis measurements of solar wind 
oxygen as thos e most likely to sample the pre-injection solar nebula. Since preliminary results 



from Genesis (IMcKeegan et al.ll2009l ; ?) (McKeegan et al. 2009, 2010) are revealing the Sun 



to be O-rich, GM07 would rule out injection of Al and Ca from a single supernova. 

In this paper, we attempted to reproduce the calculations of GM07, to apply their 
method of using oxygen isotopes to test the supernvoa injection hypothesis. We made 
necessary corrections to their method, mostly in regard to the presumed oxygen isotopic 
composition of the (post-injection) solar nebula. GM07 assumed this was identical to SMOW, 
meaning the pre-injection solar nebula had to be more 17 0-rich than almost any known 
inclusions. We presented considerable evidence that the post-injection composition was 
in fact much more 16 0-rich than that, closer to (— 60 %o, — 60 %o)- We carried out stellar 
nucleosynthesis calculations, to calculate the isotopic yields of 26 Al, 41 Ca and oxygen isotopes 
in a variety of supernova explosion scenarios, including the ID (spherically symmetric) cases 
as well as 3D (asymmetric) explosions. Because 26 Al and 41 Ca are observed to be correlated 



( ISahijpal fc Goswamil Il998l ). we also simultaneously considered injection of 41 Ca into the 
solar nebula. We then computed the shifts in oxygen isotopes and the final 41 Ca/ 40 Ca 
ratio in the solar nebula following injection of sufficient supernova material to produce the 
meteoritic ratio 26 A1/ 27 A1 = 5 x 10" 



5 



Our ID simulations largely confirm the results of GM07, that isotopic shifts are likely 
to be tens of permil and to make the solar nebula more 16 0-rich than before the injection. 
We found that injection of material from either the bulk or the explosive C/Ne burning and 
sub-explosive C burning regions of supernovae moved the composition of the solar nebula 
down the slope 1 line. Our 23 M progenitors led to isotopic shifts in oxygen which moved 
the composition of the solar nebula down the slope- 1 line, with the less energetic explosions 
producing larger shifts and time delays. The 23m progenitors, which were the most energetic 
of the 23 M cases and especially effective in producing 26 Al, generated shifts that amounted 
to only < 6%o, including a time delay of 0.7 Myr for 41 Ca to decay to its meteoritic value. 
This scenario, at least, is consistent with all the applied meteoritic constraints. If less 
than 100% of the oxygen penetrated the solar nebula material due to, for example, dust 
condensation, all but one of the ID cases are consistent with the evidence from the early 
solar system. 

We note that this conclusion differs from what GM07 infer for injection of bulk ejecta 
from 21 and 25 M progenitors. GM07 likewise found isotopic shifts downward along the 
slope- 1 line, but with a magnitude of 40 to 50 permil. It is worth noting that had GM07 as- 
sumed the same starting composition for the solar nebula, (— 60 %o, — 60 %o), that we do, then 
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they would have found the solar nebula oxygen isotopic composition to be (— 82 %o, — 82 %o) 
after injection of supernova material from an 25 M Q progenitor, and (— 81 %o, — 74 %o) after 
injection of supernova material from an 21 M & progenitor. Although these shifts are mod- 
erately large, they are down the slope- 1 line. As we established in §2.3} this would not have 
been incompatible with the meteoritic constraints, as some very 16 0-rich me teoritic samples 



in th is range are known, including CAIs in Isheyevo, at « (— 68 %o, — 66 %o) ( iGounelle et al. 



20091 ). and a ferroma gnesian cryptocrystalli ne chondrule in the CH chondrite Acfer 214, at 



~ (— 75%o, — 75 %o) (IKobayashi et al.ll2003l ). Subsequently the mass-independent fractiona- 
tion process would have shifted the nebula upward along the slope-1 line, erasing this isotopic 
shift and eventuall y producing the composition (5 17 0, <5 18 0) ~ (— 40 %o, — 40 %o) common 



to most CAIs (e.g. Iltoh et al.l I2004J ) . We conclude that the supernova injection hypothesis 



cannot be rejected based on ID models. 

Our investigation of other parameters suggest that it is even more difficult to be con- 
clusive about supernova injection. We have considered a small number of progenitor masses 
undergoing spherically symmetric collapse; in a few cases we varied other parameters, such as 
varying the explosion energy, or allowing for loss of a hydrogen envelope in a binary scenario, 
or allowing an asymmetrical explosion. In most of these cases the isotopic shifts in oxygen 
were large. Among the cases considered here, the final 5 17 values in the solar nebula varied 
from -379 % to +15 % , and the final <5 18 values varied from -379 % to +18000 % . As 
GM07 found, most of the cases where meteoritic abundances of 26 Al are injected lead to 
large (> 10 %o) shifts in oxygen isotopes. We also considered the yields in a 3D anisotropic 
explosion of a 23 M & progenitor, in the bulk ejecta and two 26 Al-rich zones analogous to 
those in the ID explosions. We find that a wide range of outcomes is possible, with oxygen 
isotopic shifts as large as hundreds of permil, or as low as < 3%o- The fact that 18 F can 
decay to 14 N instead of 18 at high temperatures makes the yield of 18 especially sensitive 
to the thermodynamic trajectory of the ejecta, which partially accounts for the spread in the 
18 yields. On the one hand, the wide range of possible outcomes makes it nearly impossible 
to state conclusively that all supernova injection scenarios can be ruled out. On the other 
hand, the wide range of possible outcomes seems to imply a degree of fine tuning so that the 
oxygen isotopic shifts in the solar nebula were not large, especially for the 3D case. 

We conclude that the hypothesis, that the 26 A1 in the solar nebula was due to super- 
nova material injected into the Sun's protoplanetary disk, can still be made compatible with 
meteoritic constraints, under two scenarios. The first is that the injected supernova material 
came from either the bulk ejecta, or from a region in a supernova that experienced ther- 
modynamic conditions like the subexplosive C burning zone. The latter is more physically 
likely. With a 0.7 Myr time delay, the injection would have moved the solar nebula oxygen 
isotopic composition from (<5 17 0, 5 18 0) ~ (— 60 %o, — 60 %o) to a more 16 0-rich value along 
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the slope one line. All but one of our explosions produce movement along the slope 1 line. 
We produce shifts as small as ~ (— 63 %o, — 63 %o), which would have produced an accom- 
panying meteoritic ratio 41 Ca/ 40 Ca = 1.4 x 10 -8 . Subsequent mixing of rocky material with 
a 16 0-depleted reservoir would have then moved the composition of meteoritic inclusions up 
the slope-1 line to values (5 17 0,5 18 0) ~ (— 50%o, — 50 %o), consistent with primitive CAIs, 
and further up the slope-1 line with time. 

The second scenario is one in which only dust grains are injected in to the protoplanetary 
disk, and very little of the supernova oxygen condenses into dust grains. If only the most 
refractory grains such as corundum were injected, then potentially x <C 10 -5 , and the isotopic 
shifts would be negligible 1 %o), for nearly all the cases considered here. It is worth noting 
that Ca is equally refractory to Al and is likely to condense from supernova ejecta under 
the same conditions that Al conden ses, so the meteoritic abundance of 41 Ca could still be 
matched following injection of 26 Al. lOuellette et al.l (120071 ) have calculated that only 1% of 



gas-phase ejecta are injected into a disk. If almost all of the Ca an d Al in the ejecta ar e 



locked up in large grains (radii > 0.1 fim) that are efficiently injected (lOuellette et al.ll2010l ). 
but less than a few percent of the oxygen is, then potentially all of the isotopic shifts in 
oxygen calculated here should be reduced by a factor of about 100. Essentially all of the ID 
cases considered here would then conform with the meteoritic constraints, and even some of 
the 3D cases as well. 

To summarize, we agree with GM07 that oxygen isotopes can be a powerful constraint 
on supernova injection models. Our calculations of oxygen isotopic shifts following injection 
from the bulk ejecta of ID supernovae broadly match the results of GM07. Had GM07 
assumed the same starting composition of the solar nebula that we did, and considered a 
smaller time delay between injection and isotopic closure, they would have found isotopic 
shifts for 20 — 25 M & progenitors that would not be inconsistent with meteoritic constraints. 
Our own calculations of the same case predict shifts that are similar, although smaller in 
magnitude, and which are also consistent with meteoritic constraints. The existence of an 
example that is consistent with the oxygen isotopic composition and the 41 Ca/ 40 Ca ratio of 
the solar nebula means that the supernova injection hypothesis cannot be ruled out. Because 
the nucleosynthesis of oxygen differs in asymmetric explosions, a much wider range of oxygen 
isotopes is possible in 3D explosions. Because of the contingent nature of the injection it 
becomes difficult to make any statement about the possibility that the solar nebula acquired 
26 A1 from such an asymmetric explosion. Finally, all oxygen isotopic shifts are reduced if only 
large grains are injected into the protoplanetary disk, and only a small fraction of oxygen 
condenses into large grains. Quantifying the fractionation of Al and O during injection into a 
protoplanetary disk is the focus of ongoing work by this research group. If only a few percent 
of the total oxygen is injected, then nearly all the ID explosions considered here could be 
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consistent with the meteoritic constraints on oxygen isotopes and 41 Ca abundances. We 
therefore conclude it is premature to rule out the supernova injection hypothesis based on 
oxygen isotopes. 
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Fig. 1. — : Shown is the time evolution of temperature, density, and entropy for representative 
examples of particles from the Ring and Bubble region of the 3D calculation. Each line is 
labelled with the abundance of 26 Al, 18 O, and 18 F (above each graph) in mass fraction per 
particle at the end of the simulation. 
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Fig. 2. — : Shown is a 2.0 x 10 11 cm thick slice in the x-z plane of the 3D simulation. 26 Al 
abundances are in red-tones (amount per particle in M ) and number ratio of 17 0/ 16 (top 
panel) or 18 0/ 16 (bottom panel) per particle in blue-tones. The sizes of the data points are 
arbitrarily chosen, but scale with their values. A color gradient was also used to visualize the 
different abundances per particle. The ligher colors/bigger data points correspond to higher 
abundances. Apparent is a Ring on either side of the center along the axis of symmetry, and 
further out from them the Bubbles, where the highest 26 Al abundance is found. 
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Fig. 3. — : Same as Figure [2] but with 41 Ca shown in green. The highest 41 Ca abundance is 
adjacent to the Rings, and only partially overlaps with the Ring- regions. 
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Fig. 4. — : Three isotope plot showing the shifts in the oxygen isotopes we calculate following 
injection of supernova material for different scenarios. The shifts from he 3D cases are plotted 
in green, those from the ID bulk cases are plotted in blue, the ID explosive C/Ne burning 
cases are plotted in orange, and the ID sub-explosive C burning cases are plotted in cyan. 
Indicated by the bigger black dots are SMOW at (0%o, 0%o) and our assumed pre-injection 
composition at (— 60 %o, — 60 %o)- The very large shift of the 3D bulk scenario was omitted 
for clarity. 



